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ABSTRACT 

We consider the dynamical evolution of a disk of stars orbiting a central black hole. In particular, 
we focus on the effect of the stellar mass function on the evolution of the disk, using both analytic 
arguments and numerical simulations. We apply our model to the ring of massive stars at ~ O.lpc 
from the Galactic Center, assuming that the stars formed in a cold, circular disk, and find that 
our model requires the presence of a significant population of massive (> lOOM©) stars in order to 
explain the the observed eccentricities of 0.2-0.3. Moreover, in order to limit the damping of the 
heavier stars' eccentricities, we also require fewer low-mass stars than expected from a Salpeter mass 
function, giving strong evidence for a significantly "top-heavy" mass function in the rings of stars seen 
near to the Galactic Center. We also note that the maximum possible eccentricities attainable from 
circular initial conditions at ages of < lOMyr are around 0.4-0.5, and suggest that any rings of stars 
found with higher eccentricities were probably not formed from circular disks. 

Subject headings: stellar dynamics — Galaxy: center — stars: luminosity function, mass function — 
methods: iV-body simulations 



1. INTRODUCTION 

The relative proximity of the Galactic Center (hence- 
forth GC) provides a unique opportunity for "close-up" 
study of processes that are expected to be crucial in 
the formation of galaxies and black holes, yet remain 
essentially unobservable in more distant galaxies. Re- 
cent advances in telescope technology have enabled us 
to resolve individual stars in the crowded GC environ- 
ment, and the development of adaptive optics has al- 
lowed determination of both velocities and positions of 
such stars to be made with ever- increa sing accuracy (e.g. 
Ghez et abirmfl^ IGenze] ^HII 1200^1 IGhez et a,1 .1 120051: 



the disk. INavaksliinl l)2006j) suggests that this will lead 
to a stellar mass function that is significantly more 
top-heavy (i.e. with significantly more massive stars) 
than that seen elsewhere in the Galaxy, and recent 
obser vational studies suggest that this is indeed th e 
case fNavakshin fc Sunvaevir2005HPaumard et al.ll2006j) . 
Other suggestions for the formation of these stellar rings 
exist, howe ver, most notab ly the infalling s tar cluster sce- 
nario (e.g. lGerhaTdll2001[ iMcMillan fc Porteeies ZwartI 



Paumard et al.l 120061) . These new data have presented 
several puzzles, notably the presence of a number of B- 
type stars very close (< O.Olpc) to the GC, and also 
the detection of one, and possibly two, coherent rings 
of massive O - fc B-type stars a t somewhat la rger radii 
of ~ O.lpc l|Genzel et all l2003t IPaumard et al. 2006). 
These stars are known to be young, and are therefore 
presumed to be the result of recent star formation at or 
close to the GC. The environment at the GC is vastly 
different from the typical environment of ongoing star 
formation in the solar neighborhood, with much larger 
pressures, densities and temperatures, as well as strong 
tidal forces, so study of the formation of these stars raises 
a number o f interesting issues (see also the recent review 
bvlAlexa,nder 2005). 

A popular theory for the origin of these rings of stars 
is that they formed via fragmenta tion of accretion disks 
around the central black hole (e.g.lLavi nfc Bcloborodovl 
20031 IGoodmanll200alGoodman fc Tanll2004t Navakshin 
2006(1 . The stars are assumed to form on nearly cir- 
cular orbits, as a result of gravitational instabilities in 
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l2^ iBerFkoff fc HansenI 120061) . This scenario can also 
result in a significantly top-heavy mass function, due to 
mass se gregation a nd stellar col lisions within the clus- 
ter (e.g.lBonnell fc Davics 199^: iGurkan fc Rasidl2005t 
iFreitag. Gurkan fc Rasio 20061 

In this paper we consider the dynamical evolution 
of a ring of stars such as those observed around the 
GC. Previous models of the stellar dynamics of such 
a rin g have considered a single stellar mass popula- 
tion (jNavakshin fc Cuadrall2005|), or the effect of pro- 
cesse s such as resonant relaxation l|Hopman fc Alexander! 
I2006|) . Here we investigate the effects of the stellar mass 
function on the evolution of the system, and use observa- 
tional determinations of the stellar orbital parameters to 
constrain the mass function of the stellar rings. We find 
that the distribution of eccentricities can tell us about 
both the upper and lower ends of the stellar mass func- 
tion, and discuss the consequences of this result for theo- 
ries of star formation in the GC environment. The struc- 
ture of the paper is as follows. In Section |21 we present 
a simple analytic model for the dynamical evolution of a 
mass-segregated ring of stars orbiting a massive central 
black hole. In Section Owe describe numerical simula- 
tions of such a system, and compare the results to the 
predictions of our analytic model. We then apply our 
model to the GC system (Section and derive con- 
straints on the initial conditions by comparing our model 
to recent observations. We discuss the consequences of 
our results, and the limitations of our analysis, in Section 
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|S1 and summarize our conclusions in Section |H| 

2. ANALYTIC MODEL 

Let us first consider the relaxation of a system of N^, 
stars of mass M,, orbiting in a disk around a black hole 
of mass Mbh (where Mbh ^ A^*-/V/,). The disk (or ring) 
is centered at radius Rq, with a radial width AR and a 
one-dimensional stellar velocity dispersion (7*. We make 
the simplifying assumption that the velocity dispersion 
is isotropic. This is not strictly valid, but it has been 
shown that the ratio of the radial and vertical velocity 
dispersions cannot beco me larger than 3 without the sys- 
tem be coming unstablejf ^ulsrud. Mark fc CarusdllQTlt 
iPolvac hcnko & Shukm;^ 1977j) . The relaxation time for 
such a system is given by 



™- G2M,p, In 

(e.g. iBinnev fc Tremaind Il987t iPaoaloizoii fc TerauemI 
120011^ where p* is the stellar density, In A* is the Coulomb 
logarithm and C is an order-of-unity constant that de- 
pends on the geo metry of the system. (Fo r a spherical 
system C ~ 0.34. IBinnev fc Tremain3ll987i ;i The stellar 
density is given by 



P* 



2TTRnAR X 2H 



(2) 



where H is the scale-height of the disk. We can express 
H = Ut, /fl in terms of the stellar velocity dispersion and 
the orbital angular velocity of the disk, il, and so for a 
Keplerian disk we have 



ilax ^ 
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where forb — 2tt/Q is the Keplerian orbital period at Rq, 
and we have absorbed a factor of 2 into the constant 
such that Ci = 2C. Consequently, the relaxation of the 
system is governed by 



G^iV^M^ In A* 
CiRoARto,h(jl 



(4) 



This form for da^/dt allows the velocity dispersion to 
increase indefinitely, and in reality a cooling term should 
be included. However, this cooling arises when stars at 
the high-velocity tail of the distribution begin to escape 
from the system, and is only significant once the velocity 
dispersion becomes comparable to the orbital velocity. 
In our application of this model to the GC we are only 
interested in the behavior of the system at relatively early 
times, and in this case this simpler treatment remains 
accurate. The effects of this simplification are discussed 
in Section |31 

We now extend this analysis to consider two differ- 
ent mass classes of stars. Mi and M2, with AIi > 
M2 and velocity dispersions cti and (72, stellar num- 
ber Ni and N2, and mean energies per particle Ei — 
3Mial/2 and E2 = 3M2cr|/2, respectively. Each 
mass class is subject to the "self-relaxation" described 
above, but in addition energy can be exchanged be- 
tween the classes. This type of analysis has previ- 
ously been used to study the evolution of a distri- 
bution of planetesimals around t he sun, and detailed 
calculations have been made (e.g. iStewart fc Wetherilil 
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Goldreich et alJ ^|2004^ point out that the treatment of 
this problem depends on whether the velocity dispersion 
of the light stars, 172, is greater than or less than the 
Hill velocity of the heavy stars, wh- The Hill velocity is 
defined as 

VH - r!i?H , (5) 

where the Hill radius Rn = i?o(Mi/Mbh)^/^. In the 
former case, referred to as the "dispersion-dominated" 
regime, the velocity dispersion is a good approximation 
for the speed of a single star and scattering encounters 
are well approximated by two-body dynamics. How- 
ever, in the latter case, known as the "shear dominated" 
regime, the tidal gravity of the central black hole is im- 
portant: in this regime the interactions are rather more 
subtle. 

In applying our model to the GC system (see Sections 
13 fc mi, we adopt Mbh = 3 x IO^Mq and R ^ O.lpc. 
Consequently, the Keplerian orbital speed is uk — 360km 
s~^, and the Hill velocity scales as 



= 7.3 



Ml 
25M, 



1/3 



kmt 



(6) 



Thus V}i <^ vk, so the Hill velocity corresponds to orbits 
with very small eccentricities. In our models, (72 > wh 
at all but extremely early times, so we work in the 
"dispersion-dominated" regime throughout. By consid- 
ering energy conservation in a two-body interaction, we 
can write the form for the "exchange" term as 



dEi 
"dT 



djh 
dt 



&G^NiN2MiM2 In A12 

C2RoARtoT:h^\2 



(7) 

Here ai2 = (a\ -I- (72)/2 is the mean of the two veloc- 
ity dispersions (and therefore the mean collision speed), 
lnAi2 is the appropriate Coulomb logarithm for such a 
collision, C2 is another order-of-unity constant, and the 
factor of 6 in the numerator arises because we are con- 
sidering the energy rather than the velocity dispersion. 
We can then use dE/dt — 3Mada/dt to find expressions 
for the relaxation of the two velocity dispersions: 



dai _ A^iMf InAi A^2MiM2lnAi2 cti 



dt 



llorbO'i 
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da2 _ iV2M|lnA2 Ar^AfiM2 In A12 a2 f Ei 



(8) 



dt 



A2t, 



2f-orb 



E2 



Here we have rewritten the constant terms for clarity, 
expressing them as Ai = CiRqAR/ for i — 1,2. The 
form of the exchange term can be considered as the 
product of a relaxation time (containing three powers 
of (T12, as in Equation^, a volume scaling term which 
accounts for the different thicknesses of the two disks 
((7i/tTi2), and a normalized energy difference. As men- 
tioned above, similar analyses are common in the study 
of planet formation. In this context other factors, such as 
physical collisions, are also significant, but if we neglect 
these terms we can compare the results to our model. 
The form of the solution is the same, and by comparison 

^ In the planet formation context, the "self-relaxation" is often 
referred to as "viscous stirring" . 
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to the three-dimensional analysis of lStewart fc Wetherilll 
|T98^ we see that that the constants Ci and C2 are re- 
lated by the ratio C1/C2 — 3.5. We adopt this ratio 
throughout. 

From these equations we can make some qualitative in- 
ferences about the evolution of this system. We see from 
Equation O that the "self-relaxation" of the heavy stars 
will be more rapid than that of the light stars'^, so we 
expect the transfer term to boost the velocity dispersion 
of the light stars, <72, while at the same time damping the 
velocity dispersion of the heavy stars cti. However, the 
rate at which this energy transfer occurs depends on both 
the velocity dispersions and the distribution of mass in 
the system. In Equation|Slthc damping term opposes the 
self-relaxation term when E2 < Ei, and we see that the 
evolution of the heavy stars can be dominated by the 
damping term only if N2M2 ^ A^iMi(ai2/(Ti)^. Since 
(T12 is at least of order ai, this means that the damp- 
ing term can dominate if the total mass in light stars is 
greater than the total mass in heavy stars. In the case 
of the light stars, however, the exchange term acts in 
the same sense as the self-relaxation term, so we always 
expect the light stars to relax more rapidly than they 
would in isolation. This simple two-component model 
shows that, in the absence of other factors such as tidal 
and resonant effects, the stellar mass function is the crit- 
ical factor in determining the relaxation of a stellar disk. 

In order to make detailed studies of the effect of a real 
mass function, it is necessary to extend this analysis to 
at least three mass classes. We add a third class with 
M3 < M2 , and consider three exchange terms of the form 
seen in Equation|7|(l-2, 1-3 and 2-3). In this case we see 
that the heaviest stars are damped by both lighter mass 
classes, while the velocity dispersion of the intermediate 
stars is boosted by the heaviest stars and damped by 
the lightest stars, and the lightest stars are boosted by 
both heavier mass classes. The distribution of stellar 
masses is critical to the evolution, but by considering 
three classes we can now emulate the effect of massive 
stars, which may only live for some fraction of the lifetime 
of the system. We apply this "three-class" model to the 
Galactic Center system in Section 0] 

3. NUMERICAL SIMULATIONS 

Our analytic model is rather simple, so in order to 
assess the validity of this approach we have conducted 
a number of numerical simulations. We use the N- 
body stellar dynamics code developed bv lHut fc Makinol 
l)2003(l . treating both the black hole and the stars as point 
masses. This code uses a Hermite integration scheme, 
which enables us to retain high numerical accuracy. 

Great care must be taken when using an A^-body code 
to analyze such a problem, as the energy of the sys- 
tem is overwhelmingly dominated by the central object. 
If we assume that the stars have typical separation H, 
then the ratio of the typical energy of a stellar encounter 
to that of the star-black-hole system is approximately 
(M,/Mbh)(i?/i?). Therefore, for IOMq stars orbiting a 
IO^Mq black hole in a disk with aspect ratio H/R = 0.1, 

^ Formally, this is true only if NiM^ > M|. However, this 
holds for any mass function dN/dM oc where F < 3 (the 

Salpeter slope is F = 2.35), and is therefore true for all cases con- 
sidered in this paper. 
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Fig. 1. — Evolution of the rms eccentricity for a system with 
Mbh = 3 X IO'^Mq, M, = 25M0 and A^* = 50. The three grey 
curves show the results of the Af-body simulations for three dif- 
ferent random realizations of the initial conditions, and the black 
curve is the mean of these three realizations. The heavy curve 
shows the mean calculated in the same manner, but with the simu- 
lations run with more strict error tolerances (see text). The dashed 
curve shows the best-fitting analytic model, with Ci = 2.2. 

the energy binding the stars to the black hole is some 10* 
times larger than the typical energy of stellar encounters. 
Consequently, we require very strict limits on the energy 
errors resulting from the numerical integration if we are 
to maintain accuracy in our simulations. This require- 
ment, combined with restrictions on computational time, 
limit us to modeling systems with relatively few stars, 
typically 150 or fewer. 

3.1. Single Mass Class 

We first consider a single mass of stars, in order to de- 
termine the value of the numerical constant Ci in Equa- 
tion 13 (and indeed to test whether a single constant is 
appropriate). We consider systems of 50 stars, orbiting 
around a black hole of mass 3 x IO^Mq. We set up the 
initial conditions as follows. The stars are distributed 
in a radial region between R = 0.05-0. 15pc (i.e. a ring 
with radius Rq = O.lpc and width Ai? = O.lpc), with 
uniform distributions in both radius and azimuth. The 
stars are given a Gaussian distribution in z, with scale- 
height H — Q.05R. All stars are given zero-eccentricity 
Keplerian orbital velocities in the x-y plane, with zero 
velocity in the z-direction, and this system is then inte- 
grated. 

We first consider a stellar mass of Af* = 25M0. We 
generated three sets of random initial conditions, and in- 
tegrated each for 6000 orbital periods at Rq (~ lO^yr for 
the parameters specified). Each model was computed 
using two different energy error tolerances, in order to 
check the numerical convergence. The cumulative frac- 
tional energy errors were typically 10~^ in one case, and 
10^^° in the second case. Translating the cumulative en- 
ergy error into a measure of the reliability of the results 
is not straightforward, but in both cases we consider the 
estimated energy errors to be sufficiently small. We com- 
puted the orbital elements of each particle twice every 
orbital period in order to plot the results. 

Figure n shows the evolution of the root-mean-square 
(henceforth rms) eccentricity in the simulations. The rms 
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eccentricity is evaluated from the simulations by comput- 
ing the instantaneous eccentricity of each star directly 
from the orbital elements. The discrepancy between the 
mean values obtained from two sets of simulations (dif- 
ferent error tolerances) is around 15%, and is comparable 
to the differences between different random realizations 
of the same simulation. Furthermore, this discrepancy is 
also comparable to the typical random fluctuations ex- 
pected (~ 15% for A'',, = 50), so we consider the simula- 
tions to be numerically converged. We see that the rms 
eccentricity of the stars rises from zero to around 0.15 af- 
ter lOMyr, and we are able to compare the results of the 
simulations to the analytic model presented in Section 
121 In a disk of small objects orbiting a massive central 
body, the rms eccentricity of the small bodies (stars) can 
be related to their velocity dispersion by 



VK 



(10) 



(e.g. lLissaueHll993D where vk = \/GMbh/ Ro is the Kep- 
lerian orbital speed. We performed a simple least-squares 
fit to determine the best-fitting value of the constant Ci . 
We evaluate the model fit by integrating Equation 0] nu- 
merically. The Coulomb logarithm is evaluated as the 
ratio of the maximum to minimum impact parameters. 
We assume that the maximum impact parameter is AR, 
and the minimum is 2GM^/a1 (e.g.^inncv & Trcmaine 
Il987j) . and therefore evaluate the Coulomb logarithm as 
A* = u^ARjlGM^. However, we note that the solu- 
tions do not depend strongly on the form adopted for 
A* . For the "low resolution" simulations the best- fitting 
value is C\ ~ 2.55; the "high-resolution" simulations 
give C\ ~ 1.85. However, we note that the velocity 
dispersion, and therefore the rms eccentricity, depends 
only on so this apparently large discrepancy in 

the scaling constant corresponds only to around a 10% 
uncertainty in the velocity dispersi on. By com parison, 
the more detailed analysis of lStewart fc Wetherilli (|1988t ) 
suggests that Ci lies in the range 1.9-2.8, depending on 
the degree of anisotropy in the velocity dispersion, so 
we are satisfied that our simplified analytic form for the 
"self-relaxation" term is accurate to within 10-15%. 

3.2. Two Mass Classes 

We now consider a system with stars of two different 
stellar masses, in order to test the accuracy of the ex- 
change term in our analytic model. We model the evolu- 
tion of a system with 10 x 5OM0 stars and 100 x lOM© 
stars. Thus the total mass in light, IOMq stars is double 
that in heavy, 50Mq stars, so we expect the eccentrici- 
ties of the heavy stars to be damped significantly. For 
comparison, we also consider a system with 10 x 5OM0 
only. Once again we evaluated three random realizations 
of the initial conditions for each; the typical cumulative 
fractional energy errors were ~ 10~^. However, with 
many more stars it was not practical to run a full con- 
vergence test in this case. Runs with more stringent er- 
ror constraints, which were run for much shorter times, 
suggest that these simulations have relaxed somewhat 
too rapidly, and that the accuracy of these simulations 
is around 20%. However, wc note that this comparison 
only considers the early evolution, where the velocity dis- 
persion is rising very steeply. 
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Fig. 2. — Evolution of the rms eccentricity for a system with 
Mbh = 3 X IO'^Mq, Ml = 5OM0, A^i = 10, M2 = WMq and 
N2 = 100. As in Fig. Q] the grey curves show the results of the 
Af-body simulations for three different random realizations of the 
initial conditions, and the black curve is the mean of these three 
realizations. Here the upper set of curves is for the 100 IOMq stars, 
and the lower curves are for the 10 50Mq stars. The dashed curves 
again show the best-fitting analytic model, with Ci = 2.2. 

The results of the two-class simulations are shown in 
Fig. 121 The analytic model provides an excellent fit 
to the comparison (10 x 5OM0 only) simulations and, 
as expected, the eccentricities of the heavy stars in the 
two-class model are significantly damped. The analytic 
curves are slightly below the simulated data, but we sug- 
gest that this is a result of slightly insufficient numerical 
accuracy in the simulations. Even with this caveat, the 
agreement is good both in terms of the relative eccentric- 
ities of the two mass classes and the absolute numerical 
scaling. Consequently, we are confident that our simple 
analytic model is accurate to within ~ 15%. 

4. APPLICATION TO THE GALACTIC CENTER 

We now apply our analysis to the GC system. One of 
the more surprising results of recent years was the detec- 
tion of large numbers of young stars very close to the GC, 
and these stars have now been well-studied observation- 
ally. A population of massive O and B stars (sometimes 
referred to in the literature as the GC "Hei stars") is 
known to exist in one, and possibly two, coherent ring- 
like structures at a distance o f around O.l pc from the GC 
e.g.lGenzel et al .I200.'^ Ic^z et al..,200,5:.Pa,nmard et alJ 



120061) . The young age 6MvrT lPaumardetani2006i) of 

these stars suggests that they formed at or very close to 
their current location, but the environment at the GC 
poses a significant challenge to conventional theories of 
star formation. Moreover, the detection of similar rings 
of stars in the center of M31 l|Bender at al.ll2005[) sug- 
gests that such systems may in fact be common, so their 
evolution warrants further study. 

One important goal of any theory of star formation 
is to predict the form of the initial stellar mass func- 
tion [henceforth (I)MF]. Observational limitations (pri- 
marily source confusion) limit the study of the lower-end 
of the mass function at the GC through direct observa- 
tions, but several indirect ap proaches have been ta ken 
to try to determine its form. iPaumard et al.l 1)20061) in- 
fer the form of the MF from their observed A'-band lu- 
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minosity function, and conclude that the slope of the 
MF r is in the range T = 0.85-1.35 (i.e. dN/dM oc 
M~^; in these imits the Salpeter slope is F = 2.35). 
iNavaks hin & Sunvaev (2005) argue that the integrated 
X-ray luminosity of the GC region observed by Chan- 
dra sets a limit on the total mass in young, low-mass T 
Tauri stars, of < 1O'*M0. We note, however, that young 
stars of earlier spectral type emit a much smaller fraction 
of their lu minosity in X-rays than their T Tauri counter- 
parts (e.g. lStelzer et al]l2006j) . so formally this limit only 
applies to stars l ess m assive than approximately 5Mq. 
INavakshin et alJ l)2006D consider th e mutual inte r action 
of the two stellar rings observed by I Genzel et al.l ()2003() 
and^aumard ct al. (2006), and suggest that dynamical 
interactions between the two rings would destroy the ob- 
served structure unl ess the total mass of t he rin gs was 
< IO-^Mq. Finally, INavakshin fc Cuadral (f200l argue 
that A^-body interactions within a stellar ring (without 
considering mass segregation) set a somewhat weaker up- 
per limit on the total stellar mass, ^ 3 x 1O^M0. These 
studies, combined with the knowledge th at ~ 3OOOM0 
is present in more massive e arly-type stars l|Genzel et al.l 
120031 iPaumard et alJl2006j) . are strongly suggestive of 
a significantly top-heavy MF in the GC system, with 
many fewer low-mass (< 5Mq) stars than would be ex- 
pected from a standard Salpeter MF. .Navakshin (|20Q^ 
suggested that such a top-heavy initial MF will arise nat- 
urally if stars are formed by the fragmentation of an ac- 
cretion disk around the central BH, and models of the 
"infalling cluster" scenario predict a similarly top-heavy 
MF (Giir kan fc Rasio 2005; .Freitag e t al. 20061. Here 
we investigate the formation of the GC system further 
by considering the effect of the MF on the dynamical 
evolution of a stellar ring. 

4.1. Basic Model 

In order to study the effect of the stellar MF on the 
dynamical evolution of the GC system we use the "three- 
class" model presented in Section |21 Paumar d et al.l 
l)2006|) identify 53 OB and Wolf-Rayet (WR) stars in 
the clockwise stellar ring, and we use this observation 
to fix the central mass class in our model. The OB and 
WR stars have masses in the range ~ 2O-3OM0, so we 
adopt M2 = 25Mo and N2 = 50. We model the ef- 
fect of a varying mass function by allowing the masses 
Ml and M3, and the slope of the mass function, F, to 
be free parameters in our model. We choose A'^i and 
N3 according to the number of stars prescribed by the 
MF at exactly Mi and M3, and thus the ch o ice of the 
MF slope fixe s A^i and IPaumard et all l|2006l see 
alsolBeloborod ov et al.ll2006D find the rms eccentricity of 
the stars in the clockwise ring to be 0.2-0.3. However, 
they find that the less well-defined counter-clockwise ring 
(containing around 20 stars) has significantly larger ec- 
centricities, with an rms value of around 0.6-0.7. How- 
ever, we clearly see from Fig.^that a ring of 50 x 25M0 
stars will reach a peak rms eccentricity of only 0.15 in 
4-8Myr. Assuming that the stellar orbits were originally 
circular (as expec ted from formation in an accretion disk: 
lNavakshinll2006|) . we therefore require the presence of 
more massive stars in order to further excite the eccen- 
tricity of these stars to the level seen at the GC. 

We choose three initial parametrizations of our model. 
Firstly, we choose a simple Salpeter MF (F — 2.35), with 
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Fig. 3. — Evolution of rms eccentricity in our "three class" model. 
In all three panels the dashed, solid and dotted lines represent the 
three mass classes from heaviest to lightest respectively, and in 
each case M2 = 25Mq and N2 = 50. The upper panel shows 
the evolution for a Salpeter mass function with Mi = 125Mq and 
M3 = IMq: in this case we see that the light stars damp the 
eccentricities of the heavier stars significantly. The middle panel 
shows the same mass function slope F = 2.35, but this time with 
a low-mass cutoff at M3 = SMq: we see that the damping effect 
is lessened in this case. The lower panel shows a top-heavy mass 
function, with F = 1.35: in this case the presence of many more 
massive stars results in much larger eccentricities. 
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RMS eccentricities of 25Mq stars, with 125Mq stars removed at t-3Myr 
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Fig. 4. — "End state" rms eccentricities of the 50 X 25M0 stars 
as a function of the MF index F. Mi = 125M0 and M3 = 5M0 are 
fixed, and the three lines show the rms eccentricity of the 25M0 
stars at 3.0, 6.0 and lO.OMyr (bottom to top, respectively). We see 
that a significantly top-heavy MF is required in order to excite the 
eccentricity of the 25M0 stars to the level seen in the clockwise 
system jPaumard et al. 2006), and also that the maximum rms 
eccentricity attainable is around 0.4-0.5. 

a maximum mass Mi — 125M0 and a minimum mass 
— IMq. We then consider a Salpeter slope with a 
"low-mass cutoff", adopting M3 = 5Mq. (As the total 
mass of a Salpeter MF diverges to low mass, this "cut- 
off" results in a lower total stellar mass.) Thirdly, we 
co nsider a significa ntly flatter MF, similar to that found 
bv lPaumard et aP |2006) (F — 1.35), while maintaining 
Ml = 125Mq and M3 = 5Mq 

The results of these models are shown in Fig. 13 We 
find that a standard Salpeter MF results in significant 
damping of the eccentricities of the more massive stars. 
This is not altogether surprising, as the Salpeter mass- 
function diverges to low mass, but we see clearly from 
Fig. O that the rms eccentricity of the 25M0 stars is 
damped to < 0.1 at the age of the GC system. Even 
allowing for the uncertainties both in our model and in 
the observed data, this is significantly lower than the 
observed eccentricities, and essentially rules out a stan- 
dard IMF for the GC rings if the stars were initially on 
circular orbits. When the low end of the Salpeter MF 
is truncated at a higher mass the rms eccentricity of the 
25Mq stars suffers noticeably less damping, but (as noted 
above) the rms eccentricity is still somewhat lower than 
observed, as there are not enough massive stars present 
to excite the eccentricities further. Only when a signifi- 
cantly flatter MF is adopted, resulting in a much larger 
number of 125M0 stars, do the 25M0 stars reach eccen- 
tricities of ~ 0.2 in 5Myr, as demanded by observations. 
We also note that these low to moderate eccentricities 
mean that few, if any, stars will be ejected from the ring 
over lOMyr timescales, even at low stellar masses. 

Fig. 0] shows the effect of varying the mass function 
slope on the rms eccentricity of the 25M0 stars, with 
Ml = 125M0 and M3 = 5M0. We clearly see that a 
Salpeter MF lacks sufficient massive stars to excite the 
rms eccentricity above 0.2. However, even extremely top- 
heavy mass functions fail to excite the eccentricities to 
values much greater than ~ 0.5 at the age of the GC 
system, which suggests that the more eccentric counter- 
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Fig. 5. — Evolution of rms eccentricity when mass-loss is taken 
into account. For clarity we plot only the eccentricities of the 
25M0 ("class 2") stars. The upper panel shows the case where the 
heaviest (I25M0) stars are simply removed at t = 3Myr. The solid 
line shows the case of a Salpeter MF, truncated at Afs = 5M0; the 
dashed line shows the case of a top-heavy mass function (F = 1.35), 
also with the same low-mass cutoff. The lower panel is the same, 
but in this case the heaviest stars were subject to progressive mass- 
loss as specified in Equation II II In both cases we see that a top- 
heavy IMF is required in order to produce rms eccentricities (in 
the 25M0 stars) larger than around 0.15. 

clockwise system ijPaumard et al.l2006|) may pose a prob- 
lem for our model. We return to this issue in Section lTTSl 

4.2. Mass-loss from massive stars 

This simple treatment of the GC system is unlikely 
to be valid, however, as stars born at the upper end 
of the IMF are subject to significant mass-loss during 
their lifetimes. Moreover, it is unlikely that stars of 
grea ter than ^ 50Mf7^ w ill live for longer than 2-3Myr 
(e.g.lSchaller e t al.ll992|) . The effect of a supernova blast 
wave on the dynamics of the stella r rings is negligible 
fe.g. lWheeler. Lecar fc McKelllQySD . but the removal of 
the massive stars from the system can have a significant 
effect, as they are no longer able to excite the velocity 
dispersion of lower mass stars. We use two crude approx- 
imations to account for mass-loss and/or stellar death in 
our simple model. 

Our first approach is simply to assume that the most 
massive stars are not subject to any mass-loss, but are 
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removed from the system at the end of their lifetimes. 
We adopt a lifetime of 3Myr for these stars, and model 
this by setting Ni = from t = 3Myr onwards. Our 
second approach is to include mass-loss from the most 
massive stars. This is done simply by making the stellar 
mass Ml a time-dependent function Mi{t). The details 
of mass loss from very massive stars are not well un- 
derstood, but the qualitative effect of mass-loss on the 
dynamical evolution of the system is not strongly depen- 
dent on the mass-loss rates adopted. For simplicity we 
adopt the mean mass- l oss rate for WR stars derived by 
iNugis fc Laine^ l^ml^ : 

login Ml = -5.73 + 0.881ogio Mi . (11) 

We use this form for the mass-loss rate Mi for Mi > 
25Mq, and assume that the stars' mass remains constant 
once it drops to 25M0. This results in a rather rapid 
decline in the mass of the heaviest stars, with the mass 
Mi{t) falling from 125Mo to 25Mq in 1.25Myr. We note, 
however, that this form was derived from observations of 
WR stars, the most massive of which were around SOM©, 
and so is of uncertain validity at very large stellar mass. 

The results are shown in Fig. [S] We see that the con- 
sequence of including mass-loss is, as expected, that the 
eccentricities of the observable stars are not excited as 
highly as in the case when no mass-loss was included. 
The exact eccentricities achieved are rather sensitive to 
some rather poorly-defined parameters in the model, 
namely the mass-loss rates and also the upper cutoff of 
the MF. However, in broad terms we see that the typical 
eccentricities of the observable stars are expected to be 
in the range ~ 0.1-0.3, depending on the details of the 
upper end of the MF. 

4.3. Implications for the Galactic Center system 

Taken together, these results have important con- 
sequences for the G C system. As noted above, 
iPatimard et al.1 l)2006(l find typical eccentricities in the 
clockwise ring of around 0.2-0.3, and in the counter- 
clockwise ring of ~ 0.7. Our analysis has shown that, in 
order for the observed stars to reach even the moderate 
eccentricities seen in the clockwise ring, there must be a 
significant population of much more massive (> IOOMq) 
stars present in numbers greatly exceeding those ex- 
pected from a Salpetcr MF. We also find that the total 
mass in low- mass (< SM©) stars must be less than the 
mass in OB stars in order to avoid significant damping 
of the eccentricities of the OB stars. Thus we conclude 
that the clockwise system has a rather top-heavy mass 
function, in agreement with the previous observ ations of 
INavakshin fc SunvaevI |2005) and lPaumardeTalJ l|2006() . 
and find that the observed eccentricities are in agreement 
with thos e expected from star formation in a Kep lerian 
disk fe.g. lGoodman Ta,nll2f)f)4t INava,kshinll200(Tt) . 

As noted in Section 14.11 however, the much 
larger eccentricities i n the counter-clockwise system 
UPa'umard et u]\\2(M iBd^orodov et a,].ll2003) pose a 
problem for our model. Even in a "best-case" model, 
with no mass-loss from the heaviest stars, and a very 
top-heavy MF (F = 0.85) extending to very high stel- 
lar mass (I75M0), the observable 25Mq stars only reach 
an rms eccentricity of 0.5 in lOMyr, and more realistic 
models suggest that rms eccentricities greater than 0.4 



are unlikely (see Fig.01). We therefore conclude that the 
counter-clockwise system was probably not formed from 
an initially circular disk, and suggest that some other 
mechanism must be responsible. 

As noted above, our model assumes initially circu- 
lar stellar orbits, but it may be possible to produce 
larger eccentricities if the initial configuration of the sys- 
tem has significant eccentricity. Such a configuration is 
also possibl e in the infalling cluster scen ario ( Gcrhar3 
120011 'McMillan fc Porteeies ZwartI l2003j) . Within the 
disk fragmentation scenario it may be possible to gen- 
erate eccentric initial conditions if the disk is itself ec- 
centric when it fragments. If the disk is formed by some 
individual "accretion event" , su ch as the captur e of a 
molecular cloud (as suggested by lNavaksh in°2006'l , then 
an eccentric stellar disk can form if the fragmentation 
timescale is shorter than the circularization timescale of 
the disk. Models of eccentric accretion disks have shown 
that a disk can typically remain eccentric for many or- 
bital times llSver fc Clarkdll992t IOgilviell200H) . so rapid 
fragmentation of such a disk may indeed provide a mech- 
anism for generating the large observed stellar eccen- 
tricities. We note also that eccentric initial conditions 
may retain an observable signature, in the form of a sig- 
nificantly anisotropic velocity dispersion. Observations 
of the stellar disk(s) at the center of M31 suggest sig- 
nificant, coherent, eccentricities, and models of eccen- 
tric stellar disks have b een shown to fit the observed 
data well (e.g. [T remainl llQQSl IPeiris fc Tremaind 12001 
IBender at al.ll2005fl . We note that the central black hole 
in M31 is approximately 100 times more massive than 
that in the Galaxy, so an eccentric stellar disk will retain 
a preferred eccentricity vector for a much longer time in 
M31 than at the GC. However, detailed study of this 
problem, and associated issues to do with the evolution 
of initially coherent eccentric stellar disks, is beyond the 
scope of this investigation. 

Our model predicts that few, if any, stars should 
be scattered inwards from the rings at O.lpc, and 
we therefore make no attempt to explain the ori- 
gin of the so-cal led "S-stars" (e.g. Genzel et al. 200^ 
iGhez et al.ll2005j) . which orbit much closer to the GC 
(at a radius of ~ O.Olpc). It has been suggested that 
these stars may have been scattered inwards from the 
rings by an intermediate-mass (IO^-IO^'Mq) black hole 
(Hansen & Milosavlievic 200^), but our model does not 
make any prediction in this regard. We note, how- 
ever, that such an interaction would provide an addi- 
tional means of exciting eccentricity in the stellar rings, 
and therefore may provide another means of generating 
the large eccentricities observed in the counter-clockwise 
ring. 

5. LIMITATIONS 

There are several obvious limitations to our analy- 
sis. Our model assumes that the velocity dispersion is 
isotropic, and while this assumption is not unreasonable 
it is not strictl y valid. Similar analy ses applied to plan- 
etesimals (e.g. IStewart fc Wetherilll Il988') suggest that 
in equipartition the radial velocity dispersion is roughly 
double the vertical dispersion. Our simulations are com- 
plicated by the fact that our initial conditions have zero 
velocity dispersion in the radial direction (Keplerian or- 
bits), but non-zero dispersion in the vertical direction 
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due to the finite thickness of the ring. The system does 
not reach equipartition over the timescales considered, 
but the general trend seems to be in broad agreement 
with previous analyses. This suggests that our relation 
between the velocity dispersion and the disk thickness 
is not exact, but we note that the manner in which we 
fit the scaling constant Ci enables us to account for this. 
Furthermore, we find good agreement between the model 
and the simulations, so we do not consider anisotropy in 
the system to be a significant problem. 

As mentioned in Section[51 our analytic model does not 
allow for cooling, and consequently does not permit an 
equilibrium solution. As long as the velocity dispersion 
remains small relative to the orbital speed (i.e. the rms 
eccentricity remains small) this approximation is valid, 
and this is supported by the favorable comparison be- 
tween the analytic model and the A'^-body simulations. 
This approach may result in an over-estimate of a (and 
therefore the eccentricity) if a becomes large. However, 
the largest eccentricities attained by any of our mod- 
els are ~ 0.5, so we do not consider this to be a sig- 
nificant problem. We also neglect cooling via physical 
stellar collisions, but note that at the stellar densities 
considered here this is unlikely to be significant. Sim- 
ilarly, although we consider only the "dispersion dom- 
inated" regime for our two- and three-class models, we 
note that our initial conditions (with circular orbits) lie in 
the "shear-dominated" regime. However, the relaxation 
of the system to the dispersion-dominated regime is very 
rapid, typically occurring within 100 orbits, so we are 
satisfied that this simplification does not affect the re- 
sults significantly. 

As noted in Section 13.21 the convergence of our "two- 
class" numerical simulations is rather marginal, with 
tests indicating that the relaxation of the system is some- 
what too fast. Unfortunately it is not practical to run 
more stringent simulations with this type (Hermite) of 
numerical algorithm, and is it not clear if more efficient 
algorithms (such as tree-codes) will provide sufficient ac- 
curacy to study this problem. Our tests suggest that 
the accuracy of our simulations is likely no better than 
±15%, but we are satisfied that this level of uncertainly 
does not affect the qualitative results of the model. We 
also note, as above, that our treatment of mass-loss from 
massive stars is somewhat arbitrary. However, the two 
models chosen span a significant fraction of the available 
parameter space, and more extreme parametrizations do 
not alter our conclusions significantly. 

As noted in Section 12 our model neglects other 
dynamical effects such as tidal and resonant effects, 
and also treats the black hole as a point mass with 
a Newtonian potential. At the large radius consid- 
ered here (O.lpc) the timescale for relativistic preces- 
sion of orbits is much longer than the age of the sys- 
tem, so we can safely neglect relativistic effects. (The 
expected precession rate is a few degrees in lOMyr: 



iWeinberg. Milosavlievic &: Ghed l2005i ) Moreover, the 
non-Keplerian compo nent of the potential due to eit her 
remnant black holes l)Miralda-Escude fc Gouldll200(]D or 
other stellar population s at the GC (e.g. iGenzel et alJ 
120031 IGhez et all 1^005") is not expected to be signifi- 
cant. The effect of reso nances on the GC syst em has 
recently been studied bv iHopman fc Alexander! (|2006), 
who find that resonant relaxation can in fact dominate 
over the type of uncorrelated two-body interactions con- 
sidered here. They note, however, that the stellar disks 
observed at the GC are sufficiently young as to be un- 
affected by resonant effects. The fact that our model, 
which neglects resonant effects, provides a good fit to 
our numerical simulations supports this conclusion. It 
may well be that resonant effects will dominate the fu- 
ture evolution of the stellar rings at the GC, but in the 
early stages of evolution considered here they do not. 

6. SUMMARY 

In this paper we have considered the dynamical evo- 
lution of rings of stars around a massive black hole. 
Through analytic arguments and numerical simulations 
we have constructed a model for the evolution of a disk 
of stars of different masses, and shown that the stellar 
mass function is the dominant factor in determining the 
evolution of such a system. We have then applied our 
analysis to rings of stars observed to orbit the Galac- 
tic Center system. We find, in agreement with previous 
studies, that the total mass in low-mass (< 5Mq) stars 
must be significantly lower than expected from a Salpeter 
mass function, and also find that a significant population 
of massive (> IOOMq) stars must have been present in 
order to pro duce eccentricit i es in t he range 0.2-0.3, as 
observed bv iPaumard et all l)2006|) . However, we find 
that dynamical relaxation alone is unlikely to produce 
rms eccentricities larger than ~ 0.4 in the GC system. 
Consequently we conclude that rings with larger eccen- 
tri cities, such as the cou nter-clockwise system observed 
bv lPaumard et al] lj2006j) , are unlikely to have originated 
in a circular disk, and suggest that some other dynamical 
process must be responsible for such systems. Alterna- 
tively, we suggest that star formation by fragmentation 
of an eccentric accretion disk could produce the observed 
eccentricities. 
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